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Abstract 

In this paper we numerically demonstrate that massless two-dimensional QCD 
is not integrable. To this aim, we explicitly solve the 't Hooft integral equation for 
bound states by an adaptive spline procedure, and compute the decay amplitudes. 
These amplitudes significantly differ from zero except in all cases in which the decay 
also produces a pion. 
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1 Introduction and Conclusions 



Integrability imposes very strong constraints on a quantum field theory. The scattering 
amplitudes become purely elastic and pair production is suppressed. These features are 
rather unrealistic from the point of view of four- dimensional quantum field theories of 
elementary particles; and indeed, integrability only exists in two dimensions. There it is 
however ubiquitous, and allows in many cases the determination of the exact S-matrix. 

One of the most interesting two-dimensional model quantum field theories is QCD2. 
Pure Yang- Mills theory in two dimensions is essentially trivial due to the absence of trans- 
verse degrees of freedom; however, adding minimally coupled fermions results in a model 
which is believed to display some of the realistic features of the four-dimensional theory. 
Indeed, as 't Hooft |] first showed in the leading approximation (N being the order of 
the gauge group), the model contains a tower of mesons built from quark-antiquark pairs, 
with nearly linearly rising meson masses. No free quarks are present. This furnishes a 
simple yet explicit example for confinement. 

Recently, it was argued in a series of papers that the model should be integrable in the 
case of massless bare quarks || for any N. This conclusion was reached from the study of a 
bosonized version of the model, in which seemingly an infinite number of conserved charges 
existed (cf [p] for a review of this and related work). An ansatz for an exact, factorized 
S-matrix was proposed in [|J. 

The assertions contained in the above work are very strong, and clearly warrant verifi- 
cation in a more direct, and controlled fashion. This may be accomplished by going back 
to 't Hooft's large N limit. Indeed, the large N limit has been a very useful tool to test the 
S-matrix of integrable vector- like theories like the O(N) model or the U(N) Gross- Neveu 
model ||. The crucial idea is to investigate the occurrence of pair production ||. 

In the following we carefully compute the decay amplitudes in the chiral 't Hooft model. 
This requires a numerical solution of 't Hooft's bound state integral equation for which no 
general analytic solution is known to date. We present in this paper an adaptive spline 
procedure, which allows to compute numerically exact solutions (both in the massless 
and the massive case) with negligible effort. The numerical code is made available to 
the interested reader; it outperforms the algorithms which were previously applied to this 
problem. The method most frequently used in the past is based on the use of Chebyshev 
polynomials and was developed in J7|. 

We definitely find non-zero probabilities for pair production to first order in 4*. This 
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result renders inconceivable the integrability of chiral two-dimensional QCD with gauge 
group SU(N). 

In contrast, we establish the complete on-shell decoupling of the massless N = oo 
pion. While expected on general grounds, the result displays a mathematically intriguing 
property of 't Hooft's bound state problem. The decoupling is satisfied to a very high 
numerical precision, which illustrates the quality of our method. 



2 Decay Amplitudes in the 't Hooft Model 



QCD 2 with gauge group SU (N) is described by the Lagrangian 

C = -^Tr F^F^ + Jii-fDn - m)V, (1) 

where the field strength is F^ = d fl A u — d u A^ + A v ] and the covariant derivative is 
Dfj, — dfj, + iA^. The gauge potential A^ is given by N x N hermitian traceless matrices. 
In the light-cone gauge and the limit iV — > oo, 't Hooft has derived the Bethe-Salpeter 
equation for meson bound states [I]: 

f m = ^(x) - £ dyj^y 2 m- (2) 

Here <f>(x) is a light-cone wave function with meson mass-square eigenvalue /i 2 and bare 
quark mass-square 7. In the chiral limit one has 7 = 0. V denotes principal value integra- 
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tion and all mass-squares are in units of — . 



At N = 00 the mesons are trivially stable since decay processes are suppressed by at 
least a factor of -k. One therefore needs the amplitude to order for a meson decaying 
into two less massive mesons. This was first worked out in |§; the amplitude A(uj) is 

1 X X — UJ 1 /" X x — UJ 

A(i,fi,f 2 ;uj) = / dx4>i{x)(j) h {-)^ h {- ) - - / dx</>i(x)$ fl (-)<i> f2 {- ). 

1 — UJ JO UJ 1 — UJ UJ J uj UJ 1 — UJ 

(3) 

Here fa, 0/ 2 are the wave functions of the initial and first and second final mesons, 
respectively. The quark and the antiquark of the initial meson go to the second and first 
final meson, respectively. The quark-antiquark-meson vertex function <&(x) (with x ^ [0, 1]) 
is related to the wave function (j>(x) through 
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w ± = ?T72 ' ( 5 ) 



and the on-shell values of the kinematic parameter uo are 

l4 + n\ ~ H% =F + -^/J 2 ~ 4 /^/4 

"2^1 

where cj + and c<j_ correspond to a right moving and left moving final state f\. Clearly pair 
production is only kinematically possible if /x, > /x^ + /x/ 2 . In the special case of identical 
final states f\ = f% = f, eq. (0) simplifies to 



- T -\ 1- A. (6) 

Note that these are partial decay amplitudes, which contain all the dynamical infor- 
mation. To get the full amplitudes entering the decay probabilities we have to add 
the partial ones. Thus we add the amplitudes for the quark of the initial particle go- 
ing to the first or second final state. One then finds for the total on-shell amplitude 
A = (1- {-\Y^h+-h){A{iJ l J 2 -uJ + ) +A{iJx,fcu-)). Here a = +1 or a = -1 for a 
state with even or odd parity. Therefore half of the decays are forbidden due to the parity 
symmetry. Below we will only consider allowed decays. 



3 The Numerical Method 

The integral equation eq. (fj) is too singular to be solved by direct discretization. One is 
practically forced to introduce a set of smooth test functions, for which the principal value 
integral in eq. (||) can be computed exactly. Besides being smooth, an ideal test function 
has to be sufficiently flexible in order not to impose its own analytical structure onto the 
solution of the equation. Cubic splines (which are also much used in computer graphics) 
provide an excellent approximation scheme for difficult eigenvalue problems as the present 
one. The power of this approach is very often overlooked. 

Our idea is to search the best solution of the 't Hooft equation among the spline 
functions 0(x) on the interval [0, 1]. The interpolating cubic spline is uniquely defined by 
its function values 0j = 4>{ x i) on a se t °f node points x i: i — 1,N with Xi < x i+ i (and by 
the derivatives at x = and at x = 1, c/0). Standard algorithms p[] ,[]T0| allow to compute 
the function which, on the interval [x{, is a third-order polynomial and which passes 
smoothly from any interval into the adjacent ones (intervals are patched together in such 
a way that </> and its first two derivatives are continuous at x^. 

For a given set of parameters (the values of fa, i = 1, N, the derivatives dip / dx\ x= o and 
d(f)/dx\ x= \, as well as /x 2 ), the eq. (0) will not be satisfied, but we can choose the parameters 
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such that the discrepancy between the l.h.s. and the r.h.s. of eq. (fj) is minimized. Using 
a standard conjugate gradient algorithm for minimization in high dimensions, we easily 
find different sets of parameters which let this discrepancy disappear for all (numerical) 
intents and purposes, fi 2 is one of the parameters, the minimization thus effectively solves 
the eigenvalue problem. 

As presented, the spline approximation imposes a locally polynomial structure, which 
is inappropriate especially close to x = and x — 1, and particularly so in the massive case 
(7 7^ 0), where the wave functions are known to behave like <p(x) ~ x 13 and 4>(x) ~ (1 — x) 13 
with 7T/9 cot 7T/3 =1 — 7. It may seem, therefore, that very large values of N have to 
be used in order to get a satisfactory solution. To avoid this unwanted feature, we also 
include the values of the Xi in our set of parameters, which the minimization algorithm 
strives to optimize. By doing this, we allow the program to choose large interval sizes where 
a polynomial approximation is locally appropriate, and small sizes elsewhere. We usually 
increase the precision of the approximation by one order to several orders of magnitudes at 
fixed N. 4 to 5 significant digits of the eigenvalues may be obtained with as few as iV = 10 
nodes. 
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Figure 1: Superposed first and 13th excited state of the massless 't Hooft equation. 
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In fig. (1) we give examples of the first and the 13th excited state of the massless 
equation (0) (all the other states may be obtained with the help of the available program). 
Note that the ground state is given by the analytically exact solution = 1. We observe 
numerically that the higher excited states states are very well described by cosines, with 
the exception of the endpoints x = and x = 1 and a slight x-dependent phase-shift when 
moving away from the midpoint x = \. 

4 Non-Integrability 

We have computed the amplitudes for all decays involving states up to the 14th excited 
state and repeated this procedure for the massive case. The on-shell amplitudes, while 
clearly non-zero in general, display a number of features that would be interesting to 
explain. Some of these features differ significantly between the massless and the massive 
case. For example, the amplitude for the decay (i, fi, / 2 ) involving a highly excited state 
i » 1 going to two low- lying states f± and f 2 goes - for increasing i - rapidly to zero in 
the massless case while it reaches a maximum in the massive case. 

Examples of our calculations can be found in fig. (2) where we display the decay am- 
plitude A(i, fi, $2'-, oj) vs lo, for different initial and final states. In figs. (2a) and (2b) we 
consider symmetric decays (i, f, /) in which the on-shell parameters u± satisfy lu + +u_ = 1, 
and where the two decay amplitudes are trivially identical. 

Note that all these decays are parity-allowed and that we obtain A(u> + ) = A{ujS). We 
have carefully checked the numerics, and clearly obtain non-zero decay amplitudes proving 
our claim: QCD 2 is not integrable. 

It should also be mentioned that the first few possible decay amplitudes are rather small 
(see e.g. (5-1-1)). Fig. (2b) (for the decay (13-3-3)) shows that some decays nevertheless 
have big amplitudes. This seems to be the case whenever we are close to threshold (i.e. when 
the momentum of the decay products is small in the rest frame of the decaying particle. 
Finally, we show in fig. (2c) a case of an asymmetric decay (9-3-1), in which eq. (|5p gives 
two non-trivially related numbers for the on-shell value of us. Nevertheless, we find that 
the decay probabilities are equal to very high precision, illustrating again the quality of 
the numerical solution. Let us mention that the same features are found in the massive 
case, which is not further discussed here. For the parity-forbidden decays we find non-zero 
probabilities of equal magnitude, but opposite sign (for any value of the mass parameter 
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7): 

A(i, A, h\ "+) = (-l) 1+< ^ + ^(i, /i, / 2 ; (7) 

which is a convolution identity that is non-trivial except for identical final states fi = f% 
(where u> + = 1 — ujJ). It is expected on physical grounds: The probability for fi going to 
the right or the left has to be the same. 



5 Pion Decoupling 

An intriguing by-product of the present work is the "experimental" fact that the pions 
cannot be produced through the decay of massive states: A = as soon as fi or / 2 are 
pions. The existence of massless pions in the 't Hooft model is a consequence of chiral U (1) 
symmetry as shown in ||11|| . Although in ([I]) both SU(N) and U(l) vector symmetry are 
gauged, the large N limit takes the U(l) gauge coupling to zero in a way that suppresses 
the chiral U(l) anomaly. The simplest argument for this is that the anomaly receives 
contributions from Feynman diagrams with one fermion loop, which by standard large N 
counting is 0(1/N). Conservation of the vector U(l) current in two dimensions allows 
writing J M = e^d^, and conservation of the chiral U(l) current = e^J 1 * then implies 



the existence of a free massless boson |12| . Thus the pion can be considered an N = oo 
artifact, and the above argument gives an indirect explanation for the observed on-shell 
decoupling. 

However, it is not at all clear in which way the decoupling manifests itself directly in 
the mathematical description of the 't Hooft model. The pion wave function (f) n is the only 
known exact solution of the chiral (7 = 0) bound state equation (Q): 

4>n(x) = 1 with jit = and $,r(x) = — -. (8) 

x[l — X) 

If, say, the second final state is a pion (f 2 = 7r) we thus expect the amplitude © to be 
zero for any initial state i and first final state f\ with /i« > jif x : 

A{u + ) = / dx )[ + -] dx U^hi—) = 0, (9) 

JO UJ + UJ + — X X — 1 UJ + Juj + LU+ 

if we choose the on-shell value (||): 

u + = 4- (10) 

We indeed verified the validity of eq. (|||) for all possible decays of the first 14 excited 
states, using our numerical method. Examples are shown in fig. (3a) (for the decay (5-4- 
0)) and fig. (3b) (for the decay (13-3-0)). In the last example, the numerically computed 
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on-shell value uo differs from the zero of the function A(u) by about 10 -4 , which is a very 
good precision considering the order of the initial excited state, and the small number of 
parameters in our spline approximation (N = 60 was used). 



v V V 
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Figure 3: Amplitudes for two sample decay processes, involving pions: (a) (5-4-0) and (b) 
(13-4-0). The vertical bar marks the on-shell value uj + . It is seen that the pions decouple 
on-shell. 

We checked that the amplitudes become non-zero as soon as the quarks acquire a mass. 
It is a challenging problem to give a mathematical proof of the convolution formula (H), 
starting from the bound state equation @. 
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